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Abstract. Cosmic rays fill up the entire volume of galaxies, providing an important source 
of heating and ionisation of the interstellar medium, and may play a significant role in the 
regulation of star formation and galactic evolution. Diffuse emissions from radio to high-energy 
7-rays (> 100 MeV) arising from various interactions between cosmic rays and the interstellar 
medium, interstellar radiation field, and magnetic field, are currently the best way to trace the 
intensities and spectra of cosmic rays in the Milky Way and other galaxies. In this contribution, 
I describe our recent work to model the full spectral energy distribution of galaxies like the 
Milky Way from radio to 7-ray energies. The application to other galaxies, in particular the 
Magellanic Clouds and M31 that are detected in high-energy 7-rays by the Fermi-LAT, is also 
discussed. 

Keywords, radiation mechanisms: general, radiation mechanisms: nonthermal, (ISM:) cosmic 
rays, ISM: magnetic fields, galaxies: ISM, (galaxies:) Local Group, (galaxies:) Magellanic Clouds, 
gamma rays: observations, infrared: galaxies, radio continuum: galaxies 



1. Introduction 

The luminosity of a star-forming galaxy like the Milky Way (MW) is dominated by 
the relatively narrow frequency range of the spectral energy distribution (SED) from the 
ultraviolet (UV) to far infrared (FIR), which is due to stellar emission and dust reprocess- 
ing in the interstellar medium (ISM). Related to the birth and death of massive stars, 
cosmic rays (CRs) are pervasive throughout the ISM (see, e.g.. Strong et al. 2007 for 
a recent review). The diffuse emissions arising from various interactions between CRs 
and the ISM, interstellar radiation field (ISRF - the UV-FIR component of the galac- 
tic SED), and magnetic field span radio frequencies to high-energy 7-rays (>100 MeV), 
but at a lower level of intensity compared to the stellar and dust component. These 
broadband emissions are currently the best way to trace CR intensities and spectra 
throughout the MW and other galaxies. Gamma rays are particularly useful in this re- 
spect because this energy range gives access to the dominant hadronic component in 
CRs via the observation of 7r^-decay radiation produced by CR nuclei inelastically col- 
liding with the interstellar gas. Understanding the global energy budget of processes 
related to the injection and propagation of CRs, and how the energy is distributed across 
the electromagnetic spectrum, is essential to interpret the radio/far- infrared relation 
(|Helou et al. 1985, Murphv et al. 20061. galactic calorimetry (e.g.. IVolk 1989p . and pre- 
dictions of extragalactic backgrounds (e.g., [Thompson et al. 200?! [Murphy et al. 2008 ), 
and for many other studies. 
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2. Broadband Spectral Energy Distribution for the Milky Way 

The MW is the best studied non-AGN dominated star-forming galaxy, and the only 
galaxy that direct measurements of CR intensities and spectra are available. However, 
because of our position inside, the derivation of global properties is not straightforward 
and requires detailed models of the spatial distribution of the emission. In recent work by 
Strong et al. (2010)|using the GALPROP (e.g., [Strong et al. 2QQQ|IMoskalenko et al. 2QQ2I : 
also see http://galprop.stanford.edu) and FRaNKIE (Fast Radiative Numerical Kode for 
Interstellar Emission - Porter et al. 2008) codes, we have calculated the broadband SED 
of the MW for the first time. Emission by stars and dust is included, together with the 
diffuse emissions for different CR propagation models consistent with local CR data. 
Figure [1] (left) shows the broadband luminosity spectrum of the Galaxy, including the 
input luminosity for CRs for a 4 kpc halo using a diffusive-reacceleration CR propagation 
model. 'Stro ng et al. (2010)| considered a range of halo sizes (2-10 kpc). For this range, 
the relative decrease in the injected CR proton and helium luminosities is ~ 10%. For 
smaller halo sizes, the CRs escape quicker requiring more injected power to maintain the 
local CR spectrum. In addition, for larger halo sizes CR sources located at further dis- 
tances can contribute to the local spectrum, which is the normalisation condition, hence 
less power is required. Contrasting with the CR nuclei, the injected primary CR electron 
luminosity increases with Zh^ which is required to counter the increased inverse- Compton 
(IC) energy losses in the halo from the longer escape time. 



galdef ID 54_z04LMS 






5 
log (Energy, MeV) 



Neutral pions 




Total gamma rays 


0.85% 




1.6% 



Figure 1. Left: Global CR-induced luminosity spectra of the MW for a diffusive-reacceler- 
ation CR propagation model with Zh — ^ kpc. Line styles: ISRF, including optical and in- 
frared scaled by factor 10~^ (magenta solid) and components for propagation model - Cosmic 
rays (dotted lines), protons (red), helium (blue), primary electrons (green), secondary electrons 
(cyan), secondary positrons (magenta); CR-induced diffuse emissions (solid lines), IC (green), 
bremsstrahlung (cyan), 7r°-decay (red), synchrotron (black, left side of figure), total (black, right 
side of figure). Right: Luminosity budget of the MW for the CR propagation model shown in 
the left panel. The percentage figures are shown with respect to the total injected luminosity in 
CRs, 7.9 X 10^° erg s"\ The percenta ges in brackets show the values relative to the luminosity 
of their respective lepton populations (primary electrons, secondary electrons/positrons). 



Figure [T] (right) illustrates the detailed energy budget for the propagation model shown 
in Fig. [1] (left). The energy channelled into 7-ray and other secondary production from 
the CR nuclei component is only a very small fraction of the total power injected in these 
particles. For the model shown, the total of synchrotron, IC, and bremsstrahlung lumi- 
nosities for the primary electrons accounts for over half of the total luminosity injected in 
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these particles. Note that it is the IC emission at 7-ray energies that is responsible for the 
majority of the energy losses, not synchrotron radiation as usually assumed. Including 
the contribution by secondary electron/positron production, approximately half of the 
total output in 7-rays is provided by these particles, even though they count for only 
^ 2% of the injected power. Non-AGN dominated galaxies like the MW turn out to be 
very good lepton calorimeters if all the energy- loss processes are taken into account. 

3. Gamma Rays from Local Group Galaxies 

Until recently, the MW was the only galaxy that was resolved in high-energy 7-rays. 
However, observations by the Fermi Large Area Telescope (LAT) have extended the sam- 
ple of resolved star-forming galaxies to include the Magellanic Clouds (|Abdo et al. 20101 
lAbdo et al.^ OTO). with M31 also being detected ^Abdo et al. 2010p . 

Figure [2] (left) shows the background subtracted counts map for a 20° x 20° region of 
interest (ROI) surrounding the LMC. The remaining feature is extended emission that 
is spatially confined to within the LMC boundaries, which are traced by the iso column 
density contour Nu = 1x10^^ H cm~^ of neutral hydrogen in the LMC ( Kim et al. 20051) . 
The extended 7-ray emission from the LMC can be resolved into several components. The 
brightest emission feature is located near (q^j2000 7 ^J200o) ^ (05^40"^, —69°15'), which is 
close to the massive star-forming region 30 Doradus (30 Dor) that houses the two Crab- 
like pulsars PSR J0540— 6919 and PSR J0537— 6910. Excess 7-ray emission is also seen 
toward the north and the west of 30 Dor. These bright regions are embedded within a 
more extended and diffuse glow that covers an area of approximately 5° x 5°. 




Figure 2. Left: Gaussian kernel (a = 0.2°) smoothed count map for the ROI about the LMC 
after subtraction of the background model for the energy range 200 MeV - 20 GeV and for a 
pixel size of 0.1° x 0.1°. Overlaid is the N(H i) contour of 1 x 10^^ H cm"^ of the LMC to 
indicate the extent and shape of the galaxy. The boxes show the locations of the 6 point sources 
(background blazars) that were included in the background model. Right: Differential average 
7-ray emissivity spectrum for the LM C. Data for model s 7ii (black dots) and 7i2 (red dots) 
(see text) and the local ISM emissivity (|Abdo et al. 20091 - grey data points) are shown. Model 
lines: solid, predicted total 7-ray emissivity computed in the framework of a one-zone model for 
7^2; long dashed, 7r°-decay; short dashed, bremsstrahlung; dotted, IC. See lAbdo et al. 20101 for 
details. 

Most of the gas in the LMC is found in the form of neutral hydrogen and helium with 
only about 1% of the total mass ionised. However, the distribution of H ll provides the 
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best fit to the observed 7-ray emission. It is characterised by a strong emission peak 
near 30 Dor which is attributed to the intense ionising radiation of the massive stars in 
this highly active region. Even accounting for the 30 Dor emission, the H ll distribution 
is still a significantly better tracer than the neutral gas of the residual emission. This 
contrasts with the MW, where the majority of the diffuse 7-ray emission appears to trace 
the neutral gas. 

lAbdo et al. 20101 tested two hypotheses for the origin of the 7-ray emission from the 
LMC: ('Hi) all 7-ray emission from the LMC is attributed to diffuse emission from CR 
interactions, and {I-L2) only emission from a disk-like component arises from CR interac- 
tions, while the 7-rays from 30 Dor originate from other sources. The emissivity spectrum 
derived using these two hypotheses for the LMC is shown in Fig. [2] (right). The integrated 
emissivity > 100 MeV for the LMC is ~ 2 — 4 times lower than the locally derived emis- 
sivity (also shown in the figure). Similar analysis for the SMC (|Abdo et al. 201Q|) and 
M31 (jAbdo et al. 2010p derive emissivities > 100 MeV that are ^ 6 - 7 and 2 times 
lower, respectively. Assuming that the proton-to-electron ratio in these galaxies above a 
few tens of GeV is similar to that in the MW and that all of the emission is from diffuse 
processes, these values point to a general picture where the intensity of the emission is 
controlled by the power injected by CR sources scaling approximately with the relative 
star-formation rates in the different systems. Future work combining the 7-ray, radio, 
and other data will result in broadband SEDs of these galaxies spanning more than 20 
decades in frequency, which can be used to investigate a variety of phenomena related to 
how energy is injected and cycled within the ISM of galaxies other than the MW. 
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